Direct Production of Tripartite Pump-Signal-Idler Entanglement in the 
Above-Threshold Optical Parametric Oscillator 
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We calculate the quantum correlations existing among the three output fields (pump, signal, and 
idler) of a triply resonant non-degenerate Optical Parametric Oscillator operating above threshold. 
By applying the standard criteria [P. van Loock and A. Furusawa, Phys. Rev. A 67, 052315 (2003)], 
we show that strong tripartite continuous-variable entanglement is present in this well-known and 
simple system. Furthermore, since the entanglement is generated directly from a nonlinear process, 
the three entangled fields can have very different frequencies, opening the way for multicolored 
quantum information networks. 
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Entanglement, which is probably the strangest of all 
quantum phenomena, is considered the most important 
resource for future quantum information technology. Re- 
cent experiments on quantum computing, storage and 
communication of information utilize different "quantum 
hardware" , such as atom clouds , quantum dots Q and 
trapped ions , all with different resonance frequencies. 
These systems will probably be used in nodes of quan- 
tum networks, implying the necessity of devising ways to 
address them without loss of quantum information. For 
networks with several nodes, multipartite entangled light 
beams will be important to carry out such tasks. 

Most of the current realizations of entangled light 
beams are implemented by combining squeezed beams 
on beam splitters 0, IE |& 0, IE 13 The beam s P lit ~ 
ter transformation is linear and does not lead to entan- 
gled beams of different frequencies. In order to produce 
multicolored entangled beams it is important to gener- 
ate them directly from a nonlinear process. In the case 
of bipartite two-color entanglement, this has been done 
very recently, in the above-threshold optical parametric 
oscillator (OPO) |H[ll[l2. 

The OPO is the best known and most widely used 
source of entangled continuous variables for quantum in- 
formation purposes Nevertheless, focus thus far has 
been on the down-converted beams it produces, usually 
overlooking quantum properties of the pump beam. Re- 
cent proposals for direct generation of tripartite entan- 
glement use so-called cascaded nonlinearities, combining 
down-conversion and sum/subtraction frequency genera- 
tion ^ij, which are not present in standard OPOs. In this 
Letter, we theoretically demonstrate that the standard 
triply resonant above-threshold OPO naturally produces 
pump-signal-idler tripartite entanglement. We show that 
the down-converted and the pump fields' noises violate 
inequalities which are sufficient for witnessing entangle- 
ment [13|. We believe this to be the simplest and most 



practical proposal of a multicolored multipartite entan- 
glement source to date. 

For tripartite systems with subsystems (k, m, n), if the 
state is partially separable, the density operator can be 
written in the form of a statistical mixture of reduced 
density operators pi^m and p^ n : 



P 



fliPi,km <o> Pi,n , 



with weights rfc > satisfying 



E 



(1) 



1. A necessary 



condition for separability of two sub-systems was demon- 
strated by Duan et al. |l5| . in the form of an inequality: 
if it is violated, there is bipartite entanglement. This 
criterion is easily checked experimentally by measuring 
second order moments of combinations of operators act- 
ing on each of the subsystems. 

The inequality presented in Ref. |l5j for the variances 
of two combinations of positions and momenta (xj,pj) 
of subsystems j = {1,2} can be readily extended to a 
combination of three subsystems [lij ]. If we define two 
commuting operators u = h\i\ + hixi + h^x^ and v = 
gxPx + .92P2 + .93P3, where the hi and <?i are arbitrary real 
parameters, for a (partially) separable state written in 
the form of eq. inequalities of the form: 



(A 2 u) + (A 2 v) > 2( 



h n g„ 



+ \hk9k + h m g m \) 



(2) 



with different permutations of the subsystems (k,m,n), 
must hold. Therefore, violations of the inequalities cor- 
responding to the three possible permutations suffice to 
demonstrate tripartite entanglement. 

For electromagnetic fields, position and momentum op- 
erators can be replaced by the field amplitude and phase 
quadrature operators, defined as functions of the cre- 
ation and annihilation operators as fij (t) = [e llpj at (t) + 

e-^ajit)] and qj(t) = i[e^aUt) - e -i ^6j(i)], where 



2 



the phase tpj of each mode is chosen from its mean value 
in order to have (g,-) = 0. In this case, p represents the 
amplitude fluctuations of the field, and q is related to 
the phase fluctuations. From the commutation relation 
[a,j,aj,] = Sjji, it follows that \pj,qj>] — In the 

present situation, we look for violations of the following 
inequalities: 



St = {A 2 (pi -p 2 )) + (A 2 ( qi +q 2 - a q a )) > 4, (3) 
S 2 = (A 2 (p + (A 2 ( qi + a 2 q 2 - go)) > 4, (4) 

S* 3 = (A 2 (p + p 2 )) + (A 2 (ai<7i + 82 - go)) > 4, (5) 

with an optimized choice of the free parameters aj, in 
order to show that all three modes are entangled, i.e. 
that the state of the full system is not even partially 
separable. 

The tripartite entangled fields are directly produced 
by a triply resonant nondegenerate optical parametric 
oscillator (OPO), composed of a \^ nonlinear crystal 
placed inside an optical cavity (a full description of field 
mean values and tuning characteristics can be found in 
Ref. |l3). The OPO is a well-known source of non- 
classical states of the electromagnetic field, both above 
and below the oscillation threshold. In this system, a 
pump photon of frequency loq is down-converted into a 
pair of signal and idler twin photons of frequencies uj\ and 
lu 2 . These fields exit the cavity and can be easily sepa- 
rated by color (pump) and polarization (signal and idler) 
in the case of type-II phase matching. Below threshold, 
signal and idler modes are in an entangled state with zero 
mean values for the electric field [l^. Above threshold, 
the parametric coupling leads to both intensity coupling 
between the three modes (this is the well-known pump 
depletion effect: a pump photon is destroyed each time a 
couple of twin and idler photons is created) and to phase 
coherence between them: the sum of the signal and idler 
field phases is locked to the pump phase as a consequence 
of energy conservation (u>i +lu 2 = uq). This leads to both 
intensity and phase correlations between the three modes 
that extend to the quantum regime, and eventually cul- 
minate in tripartite entanglement as we show below. So 
far, physicists' interest has been concentrated on the sig- 
nal and idler quantum correlations [lflj or on the pump 
squeezing [2£|. The full three- mode system has indeed 
genuine quantum properties [2l| . which are partly lost 
when one traces out the pump mode, although the signal 
and idler modes remain of course entangled [10| . 

Quantum fluctuations of the system are calculated as 
usual .22j: we start from the evolution equations of the 
operators of the three modes (So, ai,a 2 ) inside the OPO 
cavity. We write the field operators as the sum of their 
mean values and a fluctuation term and, assuming that 
the fluctuations are small compared to the mean fields, 
which is true everywhere except very close to threshold, 
we linearize these equations around the classical mean 
values |l7|. One obtains in this way six linear Langevin 



equations that enable us to calculate the evolution of 
the real and imaginary parts of the intracavity fluctu- 
ations of the three fields. If we assume that the cav- 
ity transmission factor and the extra-losses are the same 
for the signal and idler modes, the evolution equations 
can be decoupled into two independent sets 22] : two 
equations for the signal and idler difference, and four 
equations coupling the sum of the signal and idler fluc- 
tuations to the pump fluctuations. Using the input- 
output relation on the coupling mirror, one obtains the 
output field fluctuations in Fourier domain, Sp(Q) — 
[5p Q (n),6q (n),6p 1 (n),6qi(n),6p 2 (n),6q 2 (n)] T , as a 
function of the input field fluctuations. This enables us 
to determine the full 6x6 three-mode covariance matrix, 
C = (Sp(Q) Sp(— ri) T ), of the pump, signal and idler out- 
put modes, and the variance of any combination of these 
modes. The full treatment is described in Ref. [23|. 

From the calculated covariances, we derive the opti- 
mized values of the parameters a% which minimize the 
quantities Si, S 2 and S3 of Eqs. J2J - © as func- 
tions of the covariance matrix elements for the output 
field, and calculate the corresponding minimum value for 
these three quantities. We take typical experimental con- 
ditions: cavity coupling mirror transmittance for pump 
To = 10% and signal and idler beams T = 2%, and ex- 
act cavity resonance for the three modes. We can now 
study the dependence of Si, S 2 , and S3 with the normal- 
ized pump power a (power normalized to the oscillation 
threshold on resonance) and with the analysis frequency 
u> (normalized to the inverse of the cavity round trip time 
r). 




FIG. 1: (color online). Optimized sum of variances, Sf m , 
for eq. 0: <r=pump power relative to threshold, u>= analysis 
frequency, in units of f/r. 

In Fig. ^ we display the minimized value of S\ . As 
can be seen, Sf 111 is smaller than 4 in all the presented 
range of parameters, which establishes the inseparability 
of the signal and idler modes. Let us stress that the re- 
sulting violation, with the optimization of the variance, 
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FIG. 2: (color online). Optimized sum of variances, S 1 ™ 111 = 
S™ 111 , for eqs. (|1J and @: cr=pump power relative to thresh- 
old, u>— analysis frequency, in units of 1/r. 

is much stronger than that observed by tracing out the 
pump mode and looking only at signal and idler modes 
under the Duan et al. criterion |23| . In the present case, 
the measurement of pump phase increases the knowledge 
that one can obtain about the idler beam phase from the 
measurement of the signal phase. Nevertheless, the state 
can still be partially separable if the other two inequal- 
ities (eqs. 01 and |3J are not violated. The interchange- 
ability of the roles of signal and idler makes evident that 
Sz = S3. The common minimized value of this quan- 
tity is shown in Fig. We observe that it is also below 
4, implying inseparability for a broad range of values of 
analysis frequency and pump power, even though the vi- 
olation is not as strong as for S™ ln (Fig. [I}. Correlations 
between the twin beams tend to be stronger than those 
between one of the twins and the pump, since the pump 
is not generated inside the cavity. 5™ ln = S™ ln is every- 
where larger than ~ 1.7, a value obtained for a ~ 1.6. 
For this value of <r, all three fields have approximately the 
same intensities, which is in general the best situation for 
observing correlations. 

Another method to characterize the amount of entan- 
glement in a system is to study the eigenvalues of its 
covariance matrix: they provide information about the 
maximum squeezing that can be obtained from the dif- 
ferent modes by unitary transformations and about the 
maximum bipartite entanglement that can be extracted 
from these modes |'2 lj . In our case, the minimum eigen- 
value is given by the variance of p\ — f>2- The next lower 
eigenvalue is related to the combination of phases in the 
form (qi + <j2 — f3qo), where (3 is a real number. Both 
values can be strongly squeezed, at the expense of excess 
noise for the variances of their conjugate variables. 

From these two smallest eigenvalues Ai,A2 
we calculate the lo gari thmic negativity En = 
max[0, -log 2 (AiA 2 )/2] [H E|. This is a computable 
measure of the degree of bipartite entanglement of a 



system, and it is especially useful owing to its immediate 
extension to entangled mixed states. We calculate here 
the difference, E N lS , between the logarithmic negativi- 
ties for the full system and for just the signal and idler 
modes, tracing out the pump. This difference is positive 
for the full range of parameters displayed in Fig. with 
maximum values obtained for low analysis frequencies 
(ll> < 0.02/r). It is clear that quantum information 
is present in all three modes and one only recovers a 
fraction of it when restricting measurements to signal 
and idler beams. 




FIG. 3: (color online). Difference between logarithmic neg- 
ativities, -Ejv lff , for the full three modes and for only signal 
and idler modes. <r= pump power relative to threshold, ui= 
analysis frequency, in units of 1/r. 

The tripartite pump-signal-idler entanglement in the 
OPO can be observed in a broad range of frequencies and 
pump power. The correlation is, as expected, stronger for 
analysis frequencies below the cavity bandwidth T/t for 
the signal and idler modes, and for pump powers close to 
threshold, although it does not depend so much on this 
last parameter. Calculations from the covariance matrix 
show that there is a small dependence of S± and S2 on the 
cavity detunings, which is important because the locking 
of the OPO is typically done with some small detuning 
for pump and down-converted modes. If we consider the 
presence of spurious losses inside the cavity, there is a 
linear increase of the value of Si with these losses, much 
in the way observed for the intensity correlation of twin 
beams emitted from the OPO. As for S2, inseparability 
no longer occurs for lower analysis frequencies, but still 
holds for a wide range of the parameters a and lu. 

In conclusion, we have demonstrated that the stan- 
dard nondegenerate optical parametric oscillator directly 
yields tripartite entangled light beams when operating 
above threshold. Above-threshold OPOs have produced 
the highest level of intensity quantum correlations to 
date |27j. Fig. shows that they can also produce a 
very low bound for the combined phases quantum fluctu- 
ations. Thus, the magnitudes of expected quantum corre- 
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lations are among the best achievable at present. The ex- 
perimental realization of this system is much simpler than 
the proposals based on combined nonlinearities [2^, es- 
pecially considering the high degree of experimental con- 
trol achieved over the OPO. We also note that the above- 
threshold OPO entanglement renders it a possible device 
for such tasks as a tripartite teleportation network 
Moreover, it allows distribution of quantum information 
through three modes of very different frequencies, a topic 
that is attracting growing attention |l(J, |23, |30| . This is 
of practical interest, since high efficiency photodetectors 
are only available in limited ranges of the electromagnetic 
spectrum. Frequency-tunable quantum information will 
also be very useful for light-matter interfaces in quantum 
networks. 
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